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ABSTRACT 



Context. Recent studies have suggested that moving groups have a dynamic or "resonant" origin. Under this hypothesis, these kinematic struc- 
tures become a powerful tool for studying the large-scale structure and dynamics of the Milky Way. 

Aims. Here we aim to characterize these structures in the U-V-age-[Fe / H] space and establish observational constraints that will allow us to 
study their origin and evolution. 

Methods. We apply multiscale techniques - wavelet denoising (WD)- to an extensive compendium of more than 24000 stars in the solar 
neighbourhood with the best available astrometric, photometric and spectroscopic data. 

Results. We confirm that the dominant structures in the U-V plane are the branches of Sirius, Coma Berenices, Hyades-Pleiades and Hercules. 
These branches are nearly equidistant in this kinematic plane and they show a negative slope. The abrupt drops in the velocity density distri- 
bution are characterized. We find a certain dependence of these kinematic structures on Galactic position with a significant change of contrast 
among substructures inside the branches. A large spread of ages is observed for all branches. The Hercules branch is detected in all subsamples 
with ages older than ~ 2 Gyr and the set of the other three branches is well established for stars > 400 Myr. The age-metallicity relation of each 
branch is examined and the relation between kinematics and metallicity is studied. 

Conclusions. Not all of these observational constraints are successfully explained by the recent models proposed for the formation of such 
kinematic structures. Simulations incorporating stellar ages and metallicities are essential for future studies. The comparison of the observed 
and simulated distributions obtained by WD will provide a physical interpretation of the existence of the branches in terms of local or large-scale 
dynamics. 

Key words. Galaxy: kinematics and dynamics - Galaxy: solar neighbourhood - Stars: kinematics - Methods: data analysis 



1. Introduction 



Afte r the pioneering work of IProctoil ( 1869b . Kaptevn ( 1905 ) 
and iLindbladl ( 11925b . Eggen established the spatial and kine- 
matic properties of several stellar streams - hereafter classic 
moving groups - formed by stars with s imilar kinematics in the 
solar neighbourhood (see lEggenl ll996 and references therein). 
Since these structures share their kinematics with certain open 
clusters, he worked exhaustively on the hypothesis that mov- 
ing groups are a result of the dispersion (resulting from both 
external and internal causes) of stellar clusters. 

The advent of Hipparcos astrometric data led to the defini- 
tive establishment of the existence of movin g groups and to the 
recognition of subs tructures within them dChereul et alJ|l 998; 
Asiainetal . 1999a). In addition, together with Dehnen ( 1998), 
these authors first attempted to evaluate the evolutionary state 
of the members of moving groups, which is necessary in order 



to establish their origin and evolution . iDehn en ( 1998) used sub- 
samples composed of stars with different spectral types to study 
the velocity distribution o f old and young stars . As no photo- 
metric data were available. lChereul et al](ll998l) could only use 
what they called "palliative" ages to show that the age distribu- 
tion of the stars in moving groups seems to be similar to that 
observed for their whole sample. More pr ecise ages, derived us - 



ing Stromgren photometry, were used bv lAsiain et al. (1999b), 
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who observed that the velocity dispersion of several substruc- 
tures within the Pleiades moving group is compat ible with 
that e xpected for the evolution of a stellar complex (Ef removl 
Il988l) . This conclusion was restricted to groups of stars with 
ages of up to about 1 Gyr, as both differential Galactic rota- 
tion and disk heating would have dispersed older groups among 
the field stars. Using an adap t ive ke rnel and wavelet transform 
(WT) analysis. Iskuljan et all d 19991) studied a sample of 4000 
Hipparcos stars and found that the distribution function in the 
U-V plane is characterized by a few branch es that are diagonal , 
parallel and roughly equidistant. Later on. lFamaev et al.l (|2005) 
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used a maximum-likelihood method based on a Bayesian ap- 
proach to divide a sample of giant stars into several kinematic 
groups. A very wide range of ages for each of these structures 
in their corresponding H-R diagram was observed. 

The first theoretical arguments in favour of a differ- 
ent dynamic orig in of moving g roups were put forward by 
Mayor ( 1972) and lKalnajsl (1199 11) . Following this latter author, 
Dehnenl dl998 ) pointed out that orbital resonances could be the 
cause of the existenc e of most moving gr oups observed in the 
solar neighbourhood. ISMjanetalJ (11999b related the origin of 
the branches to the Galactic spiral structure, or to some other 
global characteristics of the Galactic potential, combined with 
the initial velocities of the stars. Nowadays, these dynamic or 
"resonant" mechanisms are proposed to be the most plausible 
explanation for the main moving groups. Other minor kine- 
matic structures in the solar neighbourhood, such as HR1614, 
do seem to be rem nants of a dispersed star-forming event 
( De Silva et al.l2007 ) and it has been proposed tha t other struc 



tures are related to accretion events in the Galaxy (Hel mi et al 
2006) 



The work bv lDehnenl d2000l) and lFuxl(l2000ll200lb . among 
others, focused on the origin of the Hercules structure and its 
connection to Galactic bar resonances . Furthermore, by numer- 



ically integrating test particle orbits, iDe Simone et al.l (120041 



showed that stochastic spiral density wa ves can produce kine - 
matic structures sim i lar to t hose found by Skulian et ajj d 1999b . 
Ouillen & Minchevl d2005b and IChakrabartvl d2007h have re- 
cently gone one step further and looked for the relation be- 
tween moving groups and resonances of the non-axisymmetric 
component of the Galactic potential, i.e. the Galactic bar and 
spiral structure. They have shown that by using several com- 
binations of the parameters adopted to characterize this poten- 
tial they are able to reproduce structures in the kinematic plane 
similar to those actually observed. To restrict the free parame- 
ters of the Galactic potential, all this recent work requires new 
observational constraints such as a characterization of the ob- 
served kinematic structures in the U—V plane in terms of their 
evolutionary state, their chemical composition or their possible 
dependence on Galactic position. 

Nowadays, two important observational contributions pro- 
vide new material to complement the Hipparcos and Tycho 
astrometric data: i) the CORAVEL radial velocity data for 
a significant numb er of late-type stars belonging to the 
Hippar cos catalogue (INordstrom et al .120041 for dwarf stars and 
Famaey et al.l 12005 for giant stars) and ii) the uvby-/? survey 

of FGK dwarf stars, which have allowed the derivation of ages 

i — ii — ' — I 

and metallicities (Nordstrom et al. 2004). Lastly , data on OBA- 

type stars dAsiain et al.lll999atlTorra et alj |2000) and M dwarfs 
( Rei d et al.l2002tlBochanski et al.l2005b complete an extensive 
sample of more than 24000 stars ready to be used for the char- 
acterization of moving groups. It is necessary to complement 
these new data with more powerful statistical tools to inter- 
pret both the observed and simulated data in the 4-dimensional 
space of U—V-age—[Fe/H]. 

In this paper, we characterize the observed kinematic struc- 
tures in U—V-age-[Fe / H] space by applying new statistical 
techniques with the aim of establishing observational con- 
straints that will help to reveal their origin and evolution. 



Section [2] of this paper presents the observational data we use 
together with their precisions and possible biases. Section [3] 
contains a description of the statistical method of the wavelet 
denoising (WD) we apply. Section|4]characterizes in depth the 
structures in the velocity plane. Then, Sect. [5] and Sect. [6] anal- 
yse the age and metallicity distributions of these structures. 
Finally, the main outcomes of the work and perspectives for 
the future are summarized in Sect. [7] Subsequent papers will 
further this research by using these constraints in combination 
with test particle simulations with a suitable Galactic potential 
model. 

2. Compilation of available data 

The observational data is compiled from several recent cata- 
logues and is restricted to samples with precise radial veloc- 
ities and astrometric data, as these are necessary to compute 
the U, V and W heliocentric velocity components^. In order to 
study the evolutionary state of the kinematic structures, stellar 
ages and metallicities are required. For this reason, we aim to 
cover a wide range of ages as well as spectral types and lumi- 
nosity classes. Therefore, catalogues where ages are estimated 
using photometric calibrations as well as catalogues containing 
parameters closely related to age (H ff equivalent widfbQ) are 
considered. Table Q] shows the data available in each catalogue 
and information on parameter precision. Altogether our sam- 
ple contains 24190 stars and the age of over 16000 stars can be 
used. The following is a detailed description of the sample: 

- OBA-type S t ars: S ample of 4283 stars co mpiled by 



Asiai n et al.l d!999ab (2025 A-type stars) and iTorra et al,. 
(2000) (2258 O- and B-type stars) with astrometric data 
from Hipparcos and stellar physical parameters and ages 
computed from Stromgren photometry. Most of the ra- 
dial velocities came from the radial velocity survey 



of Hipparcos earl y-type stars dFehrenbach et al.l 11987 



Grenier et al. 1999) planned to complement the kinemati- 



cally unbiased CORAVEL database for late-type stars. For 
the OB stars, when both trigonometric and photometric dis- 
tances were available, the one with smaller relative error 
was used. Relative error on photometric distance ranges 
from 14 to 25%, depending on spectral type and luminos- 
ity class. In the compilation of A stars, the Hipparcos dis- 
tance was adopted when the relative error in the parallax 
was <t k /tt < 20% and the photometric distance otherwise. 
- FGK-type Dwarfs: Sample of 13257 stars selected from the 
16682 stars in the catalogue o f the Geneva-Copenhagen 
survey dNordstrom et al.1 120041) with kinematic data, ages 
and metallicities. Proper motions come from Tycho-2 
and most of their radial v elocities were ob t ained from 
CORAVEL. As detailed in INordstrom et all d2004 . the 
Hipparcos parallax was used when <j n ln < 13% and 



1 U is the velocity component (positive) towards the Galactic anti- 
centre, V is (positive) in the direction of the Galactic rotation and W 
(positive) toward the north Galactic pole. 

2 Equivalent width of the H B emission line can be used in this case 
as an indicator of the chromospheric activity and, consequently, as age 
information. 
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Table 1. Number of stars with kinematic data, age and metallicity from each catalogue and for the total sample and error 
information. 



Data 


OBA (1, 2) 


FGK (3) 


M (4, 5) 


KM giants (6) 


total 


UVW 


4283 


13257 


863 


5787 


24190 


UVWage 


3977 


11215 


863 2 





16055 


UVW [Fe/H] 





13109 








13109 


UVWage [Fe/H] 





11215 








11215 


euxw < 2kms~' 


21% 


54% 




17% 




eu, e v , e w (kms~') 


4.5, 4.4, 3.9 


2.0, 1.9, 1.6 




4.0, 3.5, 2.8 




e age < 30% 


32% 


50% 









References. ( P lAsiain et al. i 1999a); (2) Torra et al. i 2000); (3 'Nordstrom et al. 1 2004); (4) Reid et al. 1 2002); 



(5) IBochanski et aljt2003) ; (6) IFamaev etldl 120011) . 



otherwise the photometric distance was employed. Ages 
and [Fe/H] were calculated fro m uvby-/? photometry. 
The B ayesian estimation method in J0r gensen & LindegrenI 
d2005l) was used for age computatiorH 
- M-type Dwarfs: A set of 863 stars selected from the spec- 



troscopic s urveys in Boc hanski et a l. (2005) (428 stars) and 



Reid et alj d2002l) (435 stars). The former catalogue gives 



spectroscopic distances obtained from the TiOs band, radial 
velocities from cross-correlation in the 6000-740 A range 



and U SNO-B proper motions (iMonetl 120031) . iReid et al 



(2002) use mainly proper motions and trigonometric dis- 
tances fro m the Hipp a rcos c atalo gue and radial v elocit y 



data from iGizis et all d2002h and belfosse et all (Q999). 



Both surveys give the equivalent width of the H„ line, 
which is indicative of the magnetic activity in the chro- 
mosphere and is directly related to s tellar age. Following 
the criteria in iBochanski et alj d2005l) . this sample can be 
split into two subsamples: active stars (182 young stars) and 
non-active stars (68 1 older stars) depending on whether the 
width of the Ha line is > 1A or < 1A, respectively. 
- KM-type Gia nts: A group o f 578 7 giant stars from the 



catalogue by iFamaev et al.l ( I2005I) . They belong to the 
Hipparcos catalogue and have radial velocities given in 
the CORAVEL database and proper motions taken from 
Tycho-2. Their distances were comput ed using the max - 
imum likelihood method developed by iLuri et al.l ([1996), 



which is based on a Bayesian and parametric approach. 

Figure Q] shows the distance distribution of the stars from 
the different catalogues that make up our sample. As expected, 
the volume of space sampled is highly dependent on the spec- 
tral type, which will be taken into account in our subsequent 
analysis. Distance error distributions and possible biases in the 
distance adopted are discussed in the publ ications as s ociate d 



with each catalogue. As mentioned by Skulja n et al. (1999), 



large distance errors could lead to features in the kinematic 
U-V plane being artificially radially elongated (relative to the 
point (U, V) = (0, 0)). We have checked that in our sample, dis- 
tance errors are smaller than 25%, except for 4% of stars, all of 
them KM giants, which in any case exceed 40%. 



3 In the original catalogue, ages are rounded up to the nearest 0. 1 
Gyr but in this study, ages (and up per and lower l cr confidence limits) 
rounded up to 0.01 Gyr are used 



ipper and lower lcr 
( J0rgensen 2006). 



1500 



1000 - 
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OBA 

FGK 

M 

KM giants 




100 200 300 400 500 

d (pc) 

Fig. 1. Distance distributions of stars from the catalogues that 
make up the sample. 



No s ystematic bias in the pho tometric distances is expected 



(see e.g. Nordstrom et al.l l2004l) . For trigonometric distances 



( 1/tt), iBrown et al.1 (fl 997 ) showed that a symmetric error law 



for parallaxes such as a Gaussian results in a non-symmetric 
or systematic error in distances, which leads to a biased over- 
estimated dista nce distribution. This relative bias is shown to 
be * (o-jt/n) 2 (lArenou & Lurilll999l) . This expression gives a 
maximum bias of 2% for stars with trigonometric distances of 
the FGK subsample -all of them with cr n /n < 0.1 3-, which ac- 
counts for about ~75% of all the trigonometric parallaxes in our 
sample. Moreover, the bias for the whole sample is estimated to 
be always less than 6%. The maximum likelihood method used 
to derive distances for the K M giants corrects thi s and other 
systematic trends (see Fig. 8 in lFamaey et al. I l2005h . However, 
as discussed in Sect. 14.31 other kinds of biases caused by the 
a priori parametrization of the kinematic distribution function 
cannot be ruled out. 

Heliocentric velocities and their errors have been recalcu- 
lated when the necessary data were available. For the OBA- 
type stars and KM-type giant stars, most of which have 
distances > lOOpc, the velocities have been corrected for 
Galactic differential rotation with values for Oort's constants 
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photom etry (see the error d istribution in Fig. lb of Asiain et al 



Fig. 2. Age distribution for the OBA-type stars (dotted line) and 
for the FGK-type stars (solid line). The dashed line shows the 
histogram for stars with relative errors in age of less than 30%. 



of: A = 14.82 km s-'kpc -1 and B 



- 14.8ZKms 'to 
; Whitelockll997h . 



-12.37 km s 'kpc 

dFeast & ' 

At this point it is of utmost importance to ascertain the pos- 
sible kinemat ic biases of o u r sam ple due to selection effects. 
As stated by iBinnev et al. high proper motion stars 



were preferred in pre-Hipparcos radial velocity programmes, 
which leaded to kinematically severely biased samples. During 
the last two decades, specific observational programmes for ra- 
dial velocities were unde rtaken in parallel with t he Hipparcos 
mission (CORAVEL and iFehrenbach et alJll987l) to complete 
kinematic data for stars of the Hipparcos survey. Most of the 
stars in the subsamples of FGK dwarfs and KM giants belong 
to this survey. Consequently, this bias is expected to be sup- 
press ed for 79% of the stars of our whole sample. For the OB 
stars, iTorra et al. (2000) discussed in detail a small kinematic 
bias in their subsample (see Fig. 3 therein). Numerical simula- 
tions allowed these authors to demonstrate that it had negligible 
effects on their kinematic analysis. Regarding the A-type stars. 



999al and Fig. 7 of lTorra et al]l200ri . As stated bv lTorra et al 
(2000), an overestimation in age as large as 30-50% is expected 



for highly rotating stars, that are a substantial fraction of the OB 
main-sequence stars. This, however, implies an absolute bias 
of always less than 50 Myr whose effect in our age-kinematic 
an alysis undertaken in Se ct. [5] is almost negligible. Figure 15 
in Nordstrom et al . (2004) shows the distribution of upper and 
lower lcr relative errors in age for the FGK sample. Large rel- 
ative errors in age lead to the evolutionary state of late-type 
stars being largely undetermined. Furthermore, certain biases 
in the age estimation would have non-negligible consequences 
for our study: dependence on the evolutionary models used or 
biases that arise in certain regions of the HR diagra 

Using improved calibrations, Holmberg et al. (2007) re- 
computed ages and error estimates for the INordstrom et al] 
d2004l) sample. They state that the differences between 
old (used here) and new valued are insignificant and are 
much smaller than the estimated individual uncertainties. 
Furthermore, they find an agreement b etween the new ages and 
the in dependent age determi nation by iTakedal ( 2007 ) (see Fig. 
19 of iHolmberg et al.l 120071) . All these results allow us to be 
slightly more confident of the age values and their uncertain- 
ties. However, a cut-off by error for this parameter allows us 
to work with the more reliable ages. The distribution of stars 
with relative errors in their age^l of less than 30% is shown 
in Fig. [2] This condition is fulfilled by 32% and 50% of stars 
in the OBA and FGK samples, respectively. It has to be kept 
in mind that employing this cut-off results in a change in the 
content of the working sample (e.g. uncertainties in isochrone 
ages increase when less massive main sequence stars are con- 
sidere . The most significant feature of this age histogram is 
the peak around 2 Gyr. This peak res ults from a combinatio n 
of selection effects in the sample in INordstrom et al.l (12.004ft : 
apparent magnitude, spectral type, etc. However, the lack of 
stars younger than ~ 1 Gyr -a direct result of the blue cut-off 
at b — y — 0.205- is partially corrected in our sample by adding 
in the OBA sample. 

Finally, as shown in Table [T] the [Fe/H] metallic - 
ity parameter is only av ailable for the FGK-type stars 



someof the photometric programmes on which the|Asiain et al.| in INordstro m et al.l (I2004 . The metallicity calibration by 



(1999a) compilation is based favoured stars with known pre- 
Hipparcos radial velocity data. Although a slight bias is ex- 
pected for this subsample, it only represen ts about 8% of th e 
whole sample. Finally, as was evaluated by lReid et al.l (12.002b . 
their M dwarfs catalogue -a volume limited sample- shows no 
evidence of any systematic bias (see Fig. 4 therein). On the 
other hand, the early M-type stars from lBochanski et al.l (|2005) 
were photometrically selected by color and magnitude and thus 
no kinematic bias from the selection process is expected. Thus, 
all these considerations allow us to assume confidently that no 
sample selection bias affects our kinematic study in the next 
sections. 

Figure|2]shows the age distribution for the OBA- and FGK- 
type dwarf stars in our sample. The precision of this parameter 
deserves special attention. The vast majority of the OBA-type 
stars are main sequence or moderately evolved stars and, there- 
fore, a reliable estimate for their ages is available from nvby-B 



IHolmberg et al. ( 2007b provided new [Fe/H] estimates for this 
sample. Differences between old and new estimates reach val - 
ues up to ±0.2 dex. (see Fig. 7 in Holmber g et al. 1 120071) . 



4 For instance, IPont & Eyed d2004h demonstrate that the isochrone 
dating method is subject to significant biases that arise inevitably in re- 
gions of the HR diagram where the e ffects of age on the atmospheric 
parameters become small. [Haywood ( 2006) recently studied this in 
depth and showed how biases, similar to those expected in real sam- 
ples, greatly influence the physical interpretation of the Galactic age- 
metallicity distribution. 

5 They were not availa ble when the present work was completed. 

6 INordstrom etafl d2004l) provides upper and lower limits for error 
estimates. Here we consider that a star has e age < 30% if this is true 
for both limits. 

7 For a discussion of the dependence of the degree of completeness 
of the ir sample on the cut-off by error in ages see INordstrom et alj 
J2004I) . 
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However, Reid et all ( 2007b compar ed the photometric esti- 
mate used in iNordstrom et al.l ( 120041) with values obtained us- 
ing echelle (high resolution) spectroscopy and foun d a disper- 
sion o f only ~ 0. 1 dex. In addition, as shown by lHavwood 
(2006), explicit biasing of the atmospheric parameters can lead 
to structures and spurious patterns in the age-metallicity dia- 
grams. Undoubtedly, these considerations should be taken into 
account in the metallicity study of moving groups (Sect. |6]l. 

3. Multiscale methods 

The statistical methods that we use to characterize the struc- 
tures (i.e. establish their shape, size or statistical significance) 
in the U-V-age-[Fe / H] space are: 

- Wavelet transform (WT): visualizing and detecting struc- 
tures according to their different sizes or scales (Sect. [XT). 

- Wavelet denoising (WD): obtaining a smooth distribution 
function from a point distribution via a smoothing/filtering 
treatment at different scales that eliminates Poisson fluctu- 
ations (Sect. [3~2l i. 



A detailed description of these methods can be found 
1 1994: 



Starck & Biiaoui , 
dl995l) : IStarck et al.l dl99 



Lega et al. ( 1995); Murtagh et al 



Starck & Murtagh (200?]). In this 



section, we provide a brief description of the two-dimensional 
case. In later sections these methods are applied to different 
combinations of any two of the variables considered: kinematic 
parameters U and V, age and metallicity. To perform the cal- 
culations we use the MR software^ developed by CEA (Saclay, 
France) and the Nice Observatory. It consists of a set of tools 
that implements multi-scale methods for processing ID signals, 
2D images, and 3D data volumes. 

3.1. The wavelet transform 

The WT decomposes a function f(x, y) on the basis obtained by 
translation and dilation of the so-called mother wavelet, which 
is localized in both physical and frequency space. The method 
consists of applying the correlation product between the func- 
tion and the wavelet function: 



w a (x,y) = f{x,y)® i/r\-, -) 
\a al 



(1) 



where a is the scale parameter and w a (x,y) the wavelet coeffi- 
cient. By varying a, we obtain a set of images, each of which 
corresponds to the wavelet coefficients of the data at a given 
scale. The integral of the wavelet function is equal to zero. 
Therefore, the WT analyses the overdensities and underden- 
sities of the function, assigning them positive and negative co- 
efficients respectively. A constant function would produce null 
coefficients. The key point is that the WT is able to discrimi- 
nate structures as a function of scale, and thus it is well suited 
to detecting structures at one scale that are embedded within 
features at a different, larger scale. 

It is possible to define a discrete WT which allows us to 
compute a discrete set of wavelet coefficients or a scale-related 



set of "views" of the 2-D function: the a trous WT algorithm. 
In this case, the WT performs on a grid of pixels with a bin 
size of A. A set of band-pass filters that allows structures to 
be recognized at each scale is constructed from the so-called 
scaling function. The signal co(x,y) can then be decomposed 
into a set (yv\, ...,wj,cj): 



co(x,y) = 2^ Wj(x,y) + cj(x,y) 



(2) 



where c j(x, y) is a smooth version of the original signal cq(x, y) 
and shows the details of cq{x, y) at scale of size 2 J , in units of 
the grid bin, as they have been built with a smoothing filter of 
2 J x A size. Equivalently, Cj(x, y) is obtained with a smoothing 



filter of 2 J x A size. As w 



the structures detected 



at each scale j have a size that is approximately between 2 J xA 
and 2 7+1 x A. For more detail s of the a trous WT algorithm see 
e.g. lStarck & Murtaghl (120021) . 

In the case of a set of discrete points, the function f(x, y) is 
first approximated by cq(x, y), which is obtained by smoothing 
the set of points on a grid with a bin size of A. In our case this 
is obtained from simple star counts. The choice of A depends 
on observational errors or the resolution of the data. We choose 
0.5 km s" 1 , 0.01 Gyr and 0.01 dex for velocity, age and metal- 
licity, respectively. Table [2] shows the size of the structures de- 
tected on each smooth plane and on each plane of the WT. 

Here we use the a trous algorithm for the WT with the B3- 
Spline as the scaling function. This leads to a row-by-row con- 
volution with the mask (1/16, 1/4,3/8, 1/4, 1/16), followed by 
column-by-column convolution. This algorithm has several ad- 
vantages. Due to the properties of the B3-Spline, the transfor- 
mation is isotropic and thus allows the detection of features 
with no preferred directiorj^l Secondly, it enables us to restore 
the signal cq from the set of w (Eq. [2jl easily and without losing 
any information. Finally, the complexity of the computation of 
cq(x, y) of the initial grid is proportional to the number of points 
(stars) in the sample, whereas for the WT it is O(NJ), where 
is the total number of pixels in the grid and J is the number of 
scales. This results in a very low computational cos 

3.2. The denoising method 

The data co(x,y) present statistical fluctuations related to the 
fact that the sample is finite. The noise of the signal (number 
of stars) in each initial grid bin obeys Poisson statistics. These 
fluctuations are detected as structures, especially at the smaller 
scales of the WT. The aim of WD is to ascertain the significance 
of the wavelet coefficients Wj(x,y) at each scale, which leads to 
adaptative filtering. 

If a model for the noise can be assumed, the probability that 
a wavelet coefficient Wj(x,y) is significant can be estimated. 
Then in simple thresholding methods, a detection threshold is 



http://thames.cs.rhul.ac.uk/~multires/ 



9 Note that isotropy is not defined in non-metric spaces, such as 
some of the ones used here. 

10 The computing times are 4 s for the grid construction process with 
the whole sample (24190 stars) and 1 s for the WT in a typical grid of 
400 x 400 pixels and 7 = 5 with an Intel(R) Xeon(TM) processor at 
3.0GHz. 
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Table 2. Value of A and size of the structures detected at each smooth plane (cj) and at each scale of the WT (w ■•) for velocities, 
age and metallicity. 





u,v,w 






age 






[FeJH] 






(kms- 1 ) 






(Gyr) 






(dex) 






A = 


= 0.5 






A = 


0.01 






A = 


0.01 




CO 


0.5 






CO 


0.01 






Co 


0.01 






Cl 


1 


W] 


0.75 


Cl 


0.02 


W\ 


0.015 


Cl 


0.02 


Wi 


0.015 


C2 


2 


W2 


1.5 


C'2 


0.04 


W2 


0.030 


C2 


0.04 


W2 


0.030 


C3 


4 


W3 


3 


C3 


0.08 


W3 


0.060 


C3 


0.08 


W3 


0.060 


C4 


8 


W4 


6 


C4 


0.16 


W4 


0.120 


C4 


0.16 


W4 


0.120 


C5 


16 


W5 


12 


C5 


0.32 


W5 


0.240 


C'5 


0.32 


W5 


0.240 


C6 


32 


w 6 


24 


C6 


0.64 


w 6 


0.480 


C5 


0.64 


W5 


0.480 



defined for each scale and coefficients with higher probability 
of being due to noise are rejected. But as an alternative, we 
use Wiener-like filter ing in the wavelet space , or multiresolu- 
tion Wiener filtering (IStarck & Bijaouilll994l) . This allows the 
treatment of each coefficient significance as a continuous func- 
tion and thus denoising consists of weighting each coefficient 
according to its significance. This filtering method is detailed 
in Sect. EO 

The WD method allows us to obtain a smooth distribution 
function by reconstructing the transformed data after denoising 
or filtering at each scale, which is the key point of the method. 
Thus, the denoised signal ~c{x, y) is obtained by adding the first 
J denoised scales \v] (j — I, J) to the smooth version cj(x,y): 



c(x,y) = ^ w/(x,y) + cj(x,y) 



(3) 



7'=' 



The number of scales^ J which should be used in the WD 
depends on the image or signal size. In the MR software docu- 
mentation it is stated that, in theory, this could be J — log 2 (A0, 
where N is the number of pixels of the grid in its smallest di- 
rection, but it is also suggested that, in practice, it is preferable 
to use a lower value such as J = log 2 (A0 - 1 or J — log 2 (A0 - 2. 
It seems reasonable to adopt the option of increasing J until no 
change is observed in the reconstructed denoised signal ~c(x, y), 
i.e. WD up to the scale J, hereafter called J 'plateau, where all 
signal is found to be significant (noise free). This has been the 
option followed in most of the analyses of this study. Actually, 
in most cases Jpiatem is equal to or even lower than log 2 (A0 - 1 
orlog 2 (A0-2. Nevertheless, for some specific cases, especially 
when dealing with a signal with fewer counts per pixel, we find 
that J plateau is higher than the practical values proposed. Most 
probably, in these cases, WD up to J p i a teau could imply very 
few data values in the final scales and consequently, too much 
loss of signal. These cases will be clearly flagged, as we are 
aware that conclusions derived from them would demand the 
use of larger samples for definitive confirmation. In these situ- 
ations the J up to which the WD has been carried out will be 
specified, whereas in all other cases WD has been done up to 

J plateau ■ 

WD offers several advantages since the multiscale struc- 
tures that we expect to appear could be very complex. The 



method is more straightforward than other methods used in 
this field. It offers smoothing with a unique recipe and it does 
not require additional simulations for the treatment of Poisson 
fluctuations. Furthermore, it allows an automatic local filtering 
as the denoising is carried out at several scales. The analysis 
is not restricted to one specific scale or band-width because 
all scales are visualized at the same time in the final distribu- 
tion. Finally, this method is more precise than other smoothing 
methods such as Gaussian smoothing which degrades resolu- 
tion and is shown to introduce Gaussian features into the distri- 



bution dMartmez et al.l l2005) and, therefore, is not suitable for 
the detection of structures that may be far from Gaussian. 

3.2.1. The multiresolution Wiener filtering 

In multiresolution Wiener filtering dStarck & Bijaoui 1994 ). it 



is assumed that a measured wavelet coefficient Wj, at a given 
scale j and a given position, results from a noisy process, with 
a Gaussian distribution with a mathematical expectation Wj, 
and a standard deviation cr,: 



P( Wj /Wj) = 



— 



y/2no-j 



(4) 



If the noise in the data cq is Poisson noise, as in the present 
study, the Anscombe transformation can be applied to turn it 
into a stationary Gaussian noise with unitary varia nce under 
the a ssumption that the mean value of co is large (lAnscombe 
1948b . We expect the counts associated with the structures in 
our study to be large enough for this condition to hold. After the 
Anscombe transformation is performed, the probability density 
of the coefficients Wj becomes Gaussian and the multiresolu- 
tion Wiener filtering can be applied. The noise standard devia- 
tion at each scale cry resulting from a Gaussi an noise with stan- 
dard deviation equal to 1 are calculated in IStarck & Mur tagh 
(2002) (see their Table 2). Then multiresolution Wiener filter- 
ing takes into account that the set of expected coefficients Wj 
for a given scale also follows a Gaussian distribution, with a 
null mean, as the wavelet function has null integral, and a stan- 
dard deviation S /. 



Note that the parameter n in the MR software is n = J + 1 . 



P(Wj) = 



1 



(5) 
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In the algorithm, S > is locally estimated as S 2 = s 2 - a 2 ., 
*~ J j j j 

where s 2 is the variance of wj. To estimate Wj knowing Wj, 
Bayes' theorem gives: 

P(Wj)P( W j/Wj) 1 -^Pt 



P(Wj/wj) 



where aj = 



P(wj) 



(6) 



S 2 +cr 2 



and p 2 = 



S-+CT- 



Thus, the probability 



P(Wj/wj) follows a Gaussian distribution with a mean ajWj 
and a variance /3 2 . The mathematical expectation of Wj is 
ajWj and consequently the denoised coefficients are computed 
through a linear filter with a simple multiplication of the coeffi- 
cients wj = ajWj. The complexity of the WD using this filtering 
method is O(NJ), where N is the total number of pixels of the 
grid and J is the number of scalej^l 

This and other algorithms for filtering in the wavelet space 
were compared in IStarck & Biiaoui ( 1994 ). Using an image 
with artificially added noise, they conclude that Wiener-type 
filtering provided the best results. Although it is out of the 
scope of the present study to investigate exhaustively the differ- 
ent filtering methods, we have compared the results obtained 
using a thresholding method with the multiresolution Wiener 
filtering. We conclude that the latter produces smoother distri- 
butions, without additional artifacts. 

4. Structures in velocity space 

Figure[3]shows the velocity distribution in the U—V plan J^l for 
the whole sample of 24910 stars obtained by means of the WD 
method. This figure can be compared to Fig. [4] where the ve- 
locities of the individual stars are plotted as dots. Table [3] and 
Fig. [3] show the pos itions of the ce ntre s of the classic mo ving 



groups according to lDehnenl dl998l) and lEggenldl97lUl996lf 



14 



Connections and continuity can be seen betwee n these classic 



kinem atic structures. As was first pointed out by Skuli an et al 
19991) . these connections arrange classic moving groups in 



several branches whose approximate positions can be shown 
also in Fig. [5] These superstructures (the branches), structures 
(moving groups) and their substructures are characterized in 
this section. 

4.1. Classic moving groups 

The established understanding of classic moving groups imme- 
diately leads to the following points when analysing Fig. [5] 

- The Sirius, Coma Berenices, Hyades and Pleiades moving 
groups are clearly identified but they have neither a clearly 



12 The computing time for a typical grid of 400 x 400 pixels and 
J = 5 is 4 s with an Intel(R) Xeon(TM) processor at 3.0GHz. 

13 We deal only with these two components of the velocity because 
classic moving groups are less visible in the W distribution due to the 
more thorough phase-mixing of the movement perpendi cular to the 
Galactic plane dDehnerJll998l;ISeabroke & Girmorell2007l) . 

14 The list could be far more extensive; only the most signif- 
icant groups have been chosen. Similar values for their compo- 
nents have been reported by other authors (see e.g. iMontesI 1200 ll ; 
lLopez-S antiago et al. 2006). 



Table 3. Heli ocentric velociti es of the main moving group s ac- 
cording to (l) lDehnenld 1998b and (2) lEggenldl97lHl996l and 
references therein). 





(1) 




(2) 




Moving group 


U 


V 


U 


V 




(tar 1 ) 


(kms -1 ) 


(kms-') 


(kms" 1 ) 


1 Pleiades (Stream 0) 


-12 


-22 


-11.6 


-20.7 


2 Hyades (Stream I) 


-40 


-20 


-40.4 


-16.0 


3 Sirius (Stream II) 


9 


3 


14.9 


1.4 


4 Coma Berenices 


-10 


-5 






5 NGC 1901 


-25 


-10 


-26.4 


-10.4 


6HR1614 


15 


-60 


5.8 


-59.6 


7 


20 


-20 






8 


-40 


-50 






9 


-25 


-50 






10 


50 









11 


50 


-25 






12 IC 2391 






-20.8 


-15.9 


13 


-70 


-10 






14 


-70 


-50 






15 61 Cygni 






-80 


-53 


16 ( Herculis 






-30 


-50 


17 Wolf 630 






25 


-33 



defined shape nor defined limits. Furthermore, they present 
substructure. 

- The centre of the Sirius moving group is not well-defined, 
which could explain the discrepancies in its position on 
the U-V plane found in the literature. This moving group 
seems to be better described by a branch-like shape with a 
clear extension. 

- Similar arguments can be applied to the Hercules moving 
group, although it is less prominent in Fig. 0due to its low 
density. The two peaks proposed bv lDehnenl d 1 9981) (groups 
8 and 9 in Table are faintly observed but this structure 
appears to be a continuous elongated feature. 

- The Hyades and Pleiades moving groups are clearly 
grouped together, forming just one branch. The centres es- 
tablished by Eggen and Dehnen for NGC 1901 and IC 2391 
(groups 5 and 12 in Table [3) seem to be placed slightly 
above the crest of the branch (less negative V), in a low 
density region within the distribution. 

- Groups such as NGC 1901 (group 5), HR 1614 (group 6) 
or IC 2391 (group 12) are not observed, possibly due to the 
small contribution they make compared to the main struc- 
tures. 

- New structures such as that centred at (35, -20) kms -1 , 
which is weak and does not have a well-defined shape, 
may all be considered as part of the elongation of the Sirius 
or Coma Berenices structures, rather than being necessar- 
ily linked to the nearest identified groups (1 1 or 17 in this 
case). 



4.2. Kinematic branches 



Sku lian et alj (11999b detected the existence of at least three 
long, parallel and equidistant branches in the U—V plane: 
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Fig. 3. Density field in the U—V plane for the whole observational sample obtained by WD (J p i a reau - 4). 
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Fig. 4. Stellar velocities in the U—V plane for the whole obser- 
vational sample. 



the Sirius branch, the middle branch (here Coma Berenices 
branch) and the Pleiades branch (here Hyades-Pleiades 
branch). The increased extent of our sample, in terms of both 
luminosity and spectral type, allows us to characterize these 
three branches and also the Hercules branch. In this section the 
properties of these branches in the whole sample are studied 
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Fig. 5. As Fig.[3]but with the superposition of the classic mov- 
ing groups in Table [3] the approximate trace of the branches 
and the trace of the edge line (see text in Sect. l4.2t . 



and Sect. 14.31 deals with the subsamples of different spectral 
type. 

To simplify the study of the branches, a clockwise rotation 
through an angle (3 is applied to the original (U, V) components. 
In the new coordinate system (U rot , V rot ) the branches are better 
aligned with th e horizontal axis (ex cept for their slight curva- 



19991) proposed a turn of 25° 



ture). Although ISkulian et al. 
value of /J ~ 16° is more suitable for the data we present here, 
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V rol (kms" 1 ) 

Fig. 6. Density field for the V ro , component obtained integrat- 
ing the density in the U rot -V rot plane. The approximate posi- 
tions of the four branches of Hercules, Hyades-Pleiades, Coma 
Berenices and Sirius are marked with vertical lines according 
to the maximums. 



T 




U rot (kms ') 

Fig. 8. Density field for the U rot component obtained integrat- 
ing the density in the U ro t-V ro t plane for each of the four 
branches (Jpiateau = 4). 




Ui(kms-') 

Fig. 7. Density field in the direction perpendicular to the edge 
line (U±) obtained integrating the density in the U—V plane for 
the whole observational sample (Jpiateau - 4). 



especially for the three main branches. Their approximate po- 
sitions are delineated as dotted lines with this slope in Fig. [5] 

Figure [6] shows the density along the V rot component for 
the whole sample obtained integrating the density in the U rot - 
V ro , plane. The WD has been carried out up to different / (3, 
4, 5) in order to illustrate the discussion regarding the choice 
of J in Sect. 13.21 The four branches are clearly seen and their 
approximate positions are marked with vertical lines. It can be 
noticed that J has become rapidly Jpiateau and results almost do 
not change from J to J. Our data confirm the nearly equidistant 
character of the branches. The first three branches are found 
to be separated in intervals of approximately 15 kms -1 . The 
possible dependence of this interval on spectral type or age 
is discussed in Sects. 14.31 and [5] The separation between the 
Hy ades-Pleiades and the Hercules branch is ~ 30 kms 1 . 



Skulian et all (11999) noticed the existence of a fairly sharp 



'edge line' at 30° to the U axis and connecting the extremes of 



the branches at low U and high V (see their Fig. 5). This feature 
can be evaluated by computing the density field in the direction 
perpendicular to this hypothetical 'edge line' (direction of U± 
from now on). First, several directions were scanned to look 
for the sharpest distribution: a lower inclination of about 20° 
was found here. This edge line is plotted in Fig.[5]and the per- 
pendicular density distribution along U ± is presented in Fig. 
Q Undoubtedly, the asymmetry and the long tail at positive 
U± > observed would support the existence of a sharp edge 
line. However, we want to stress that: 

- The upper part of the edge line is not empty: a considerable 
number of stars occupy this region of the U—V plane (see 
Fig. Hi. 

- Equivalent importance should be given to other similar 
'edge lines' such as the sharp distribution observed in Figs. 
|4] and [3] near the Hyades-Pleiades branch towards low U 
and low V, which is shown in Fig. [6] with the drop at the 
left of this branch. 

- Each branch presents a different density distribution near 
the edge line (Fig. [3}. This is also confirmed by examin- 
ing the density characteristics along each branch, i.e. in the 
direction of the component U rot , plotted in Fig. [8] These 
density distributions are obtained from the whole U rot - 
V rot distribution considering the central positions in V rot of 
the branches according to the maximums in Fig. [6] (-55, 
-25, -10 and 5 km s" 1 for the Hercules, Hyades-Pleiades, 
Coma Berenices and Sirius branches, respectively) and a 
width of ±4 kms" 1 . It can be seen that the slope of the 
density drop at negative U ro , is different for each branch, 
being considerably steeper for the Hyades-Pleiades branch 
in part due to the gathering of stars of the small substruc- 
tures at the very beginning of the branch. The drop be- 
comes less abrupt from Coma Berenices to Sirius to the 
Hercules branch. Moreover, the Hyades-Pleiades and Coma 
Berenices branches show a slightly asymmetric distribu- 
tion, with a longer tail of stars at positive U rot . 
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The existence of all these abrupt features in the U—V plane 
definitely rules out the classic idea of a smooth velocity field 
distribution and it even questions the scenario where several 
structures are superimposed on this smooth field. 

4.3. Spectral type analysis 

Figure [9] shows the density field in the U—V plane (left) and 
the density for the V rot component (right) for the subsamples of 
OBA, FGK, M dwarfs and KM giants. As mentioned in Sect. 
[JJ some of these subsamples have been used for moving group 
studies by different authors and using different techniques. This 
is the first time that the same statistical method is applied to all 
of them. In this sense, a direct comparison among subsamples 
deserves special attention. 

In the right column plots, the WD has been carried out up 
to different J (3, 4, 5), as in Fig.[6]for the whole sample, in or- 
der to return to the considerations about the choice of J in Sect. 
13.21 The maximum value of J presented in these plots is always 
higher than or equal to Jpiateau, which is the one for which the 
result does not change when higher scales are denoised. For the 
FGK and OBA subsamples, J becomes rapidly Jpiateau an d the 
signals almost do not change from J to J. Therefore, the signif- 
icance of the structures and branches present in these distribu- 
tions is strongly supported. However, for the subsamples of M 
dwarfs and KM giants the peaks whose existence has been in- 
dependently proved by the previous robust results here are dras- 
tically diluted when denoising up to the Jpiateau- For this and to 
avoid the loss of signal (Sect. 13. 2t . we present and discuss the 
U—V distributions with 7 = 4 for all subsamples despite this 
not being J p i ate au for all of them. 

Some of the well-known characteristics according to spec- 
tral type are observed in these figures, for example the predom- 
inance of the Coma Berenices and Pleiades groups or the ab- 
sence of the Hyades moving group for the OBA stars. As some 
of these features depend on age, they will be studied in Sect. 
[5] More importantly, the three branches of Hyades-Pleiades, 
Coma Be renices and Si r ius ar e detected in all subsamples (Fig. 
|U right). ISkulian et al.1 dl999t) . with a different rotation angle 
P, suggested a separation of 15kms -1 for early-type stars but 
20kms~ 1 for late-type stars. Here, the value of 15kms _I is 
more suitable for all subsamples. In all cases, the positions of 
these three overdensities differ by less than 3-4 km s~ 1 from the 
values derived for the whole sample (vertical lines). 

The comparison between the velocity distributions of the 
KM giants and the FGK dwarfs deserves special attention since 
they have very different space volume coverage (see Fig. [JJ| 
and the number of stars in both subsam ples is comparable. A 
similar examination was carried out by ISeabroke & Gilmorel 
(2007). We agree with these authors that the kinematic struc- 
tures are well maintained in both distributions and thus, con- 
sidering their different spatial extension, moving groups can- 
not be simple remnants of star clusters. Our method allows us 
to go further in this comparison. For instance, a clear overden- 
sity appears at (U, V) ~ (-27, -22)kms~ I in the middle of 
the Hyades-Pleiades branch for the KM giants sample which 
is not present for the FGK dwarfs. Most specially, the Hercules 



branch presents significant differences both in shape and in po- 
sition of the density maximum. Whereas for the KM giants 
this branch is seen as a clear elongated structure, for the FGK 
dwarfs there is a density maximum more localized in the U— 
V plane. Moreover, Fig. [9] (right) shows an evident discrep- 
ancy between the positions of the Hercules branch located at 

Vrot 62kms _1 for the KM giants but at V ro , 55kms _1 

for the FGK dwarf stars. Two possible explanations for this 
discrepancy are: i) errors in distance estimates or biases and 
ii) real differences due to the different volume coverage or 
Galactic position. Previous discrepancies in the positio n of 
the H ercules group are seen by checking the literature. iFux 



(2001) found Hercules centred at (U, V)=(-35, -45) km s" 



us- 



ing a sample of nearby stars (d < 
of (-42, -51) km s~' was obtained by iFamaev et al 



1 00 p c), whereas a value 
" d2005l) 



for the KM giants. These values correspond to V rot = -53 
and -61 km s _1 , respectively. Also, as mentioned in Sect. 14.11 



Dehnen ( 1998) detected two peaks in the region of this branch 
(groups 8 and 9 in Table [3]l at V rot = -59 and -55kms _1 . 
Although density values in the tail of the distribution where 
Hercules is placed are marginally significant, we will pay spe- 
cial attention to this discrepancy. 

As mention in Sect. [2] no significant bias is expected for the 
FGK stars with trigonometric parallaxes (with small relative 
errors of less than 13%). We have checked that the position 
of the Hercules branch does not change at all when separately 
considering the FGK stars with trigonometric parallax or those 
with photometric distances. Consequently, this establishes the 

existence of the peak at V lot 55 km s _1 for stars nearer than 

100-150 pc. 

More importantly, we have selected stars from the subsam- 
ple of KM giants with relative errors in trigonometric paral- 
laxes cr n ln < 10%. As mentioned, biases in trigonometric dis- 
tances are negligible for this truncation of the relative error. In 
Fig. [10] we compare the density in the V rot component obtained 
calculating th e velocity using: i) the distance estimates from 
the method in Luri et al.l d 1996b -LM distances from now on- 
(solid line) and ii) the distances from trigonometric parallaxes 
(dashed line). As expected, both curves are practically identical 
using LM and trigonometric distances since high-quality infor- 
mation on trigonometric parallaxes was used in the derivation 
of the LM distances. Remarkably, we observe that the overden- 

sity of the Hercules branch appears at V ro , 55 km s _1 , as for 

the FGK dwarfs. With this cut in cr^/n, we are in fact selecting 
stars with distances < 150 pc, which proves again the existence 

of a real kinematic structure at V rot 55 km s _1 for the nearer 

stars. Note that, by contrast, we find that the velocity distri- 
butions of the stars located in shells of LM distances further 
than lOOpc centred in the Su r0 do present always the peak 

at 62kms _1 using the LM distances (figures are omitted). 

Therefore, the situation is complex. 

We want to evaluate whether the subsample of M dwarfs 
could contribute to solve this issue, despite being aware that 



15 Note that due to its source of radial velocities, this subsample is 
restri cted to the northern equatorial hemisphere (see IFamaev et alj 
120051) and therefore, stars in these shells are not isotropically dis- 
tributed around the Sun. 
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Fig. 9. Left column: Density field in the U—V plane for the the subsamples of OBA, FGK, M dwarf stars and KM giant stars 
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Fig. 10. Density field for the V rot component obtained integrat- 
ing the density in the U rot -V ror plane for J = 4 of the WD for 
the KM giants with cr^/n < 10% (756 stars) using the trigono- 
metric distance (dashed line) and the LM distance (solid line). 




Fig. 11. Left column: Density field in the U-V plane of two 
subsamples of KM giant stars situated at different galactocen- 
tric radii obtained by WD with J - A: inner subsample (R < 
R Q - 20 pc, 1812 stars) and outer subsample (R > R e + lOOpc, 
2057 stars). Right column: Density field for the V ro t compo- 
nent obtained integrating the density in the U ro t-V ro t plane with 
different values of J for these subsamples. 



the low number of stars in this case prevents us from giving 
any conclusive statement. The V rot distribution for these stars 
(Fig.|9]l shows two peaks around the Hercules branch at V rot ~ 
-50kms _1 and -65kms~'. As this subsample is composed of 
stars from two diff erent catalogues (s ee Sect. [2]), we ascertain 
that the stars from Reid et al. (2002) (very nearby stars with 
d < 25 pc) show a very broad Hercules branch but centred at 
V rot ~ -55kms _1 , in agreement with the previous results for 




Fig. 12. Left column: Density field in the U-V plane of two 
subsamples of KM giant stars situated at different azimuths ob- 
tained by WD with J — 4: subsample towards anti -rotation (rj < 
40 pc, 1903 stars) and subsample towards rotation (rj > 160pc, 
1975 stars). Right column: Density field for the V ror compo- 
nent obtained integrating the density in the U ro ,-V ror plane with 
different values of J for these subsamples. 



nearby stars. However, the stars from lBochanski et al.l (|2005) 
(stars from selected areas mostly at 6 > and with larger dis- 
tances up to ~ 150 pc) contribute to the split of this branch into 
two peaks. Again, the importance of the mean distance and spa- 
tial distribution of the subsample considered is revealed. 

We look for possible variations of the kinematic proper- 
ties of the branches with Galactic position. At present, the KM 
giant subsample is the only one with enough stars and space 
volume coverage to undertake this study. We have selected 4 
subsamples at different locations and approximately with the 
same number of stars (1800-2000 stars). Note that the cen- 
tre of KM giants spatial distribution onto the Galactic plane is 
placed 100 pc away from the Sun in the direction of Galactic 
rotation and 40 pc towards the Galactic anti-centrj^. Stars in 
the central region are not considered in order to emphasize the 
properties of the extremes. First, we divide the sample into two 
subsamples according to their different galactocentric radii, R: 
R < R e - 20 pc (inner subsample) and R > R Q + 100 pc (outer 
subsample). Secondly, we build two subsamples with different 
T] (heliocentric Cartesian coordinate towards the direction of 
the Galactic rotation): with rj < 40 pc (subsample towards anti- 
rotation with respect to the centre of the sample) and rj > 160 pc 
(subsample towards rotation). 

The results are shown in Figs. QT| and [12] The same 
color scale is used in all cases for a clearer comparison. 
From these figures we conclude that the four branches are 



16 Note also that the observational spatial restriction affects the 
distribution of the subsamples in the direction perpendicular to the 
Galactic plane. 
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present in all regions. Second, the Hyades-Pleiades branch 
is the dominant structure, except for the region towards 
Galactic rotation where Coma Berenices has the same den- 
sity. Furthermore, a significant change of contrast among sub- 
structures inside the branches is confirmed and the density 
maximum along the branches (along the U rot component) 
varies for each region. Note for instance the substructures at 
(U,V) ~ (-27, -22) km s" 1 and (U,V) ~ (-20, -10) km s" 1 
or the Pleiades moving group. Also the shape of the Hercules 
branch changes between regions and it is more significant for 
the region through anti-rotation. All these considerations sug- 
gest a real effect of Galactic position on the shape of kinematic 
structures, independently of the possible bias in the LM dis- 
tances. 

Resuming the issue of the discrepancies in the Hercules 
branch position, we observe that this branch appears at V rot ~ 
-62 km s _1 for all these four different galactocentric directions 
of the KM giants subsample, in strong contrast with the peak 

at V ro , 55kms _1 for the central nearby stars of this same 

subsample. This makes the method for LM distance derivation 
suggestive of a possible bias entangling kinematics and dis- 
tances in a complex way. Although the Bayesian parametric 
LM approach was developed to derive unbiased distances, the 
a priori adoption of a Schwarzchild ellipsoid for the velocity 
distribution functions can directly affect the distance estimate 
of a given kinematic group or branch. Further work along these 
lines requires the revision of the LM distances and the applica- 
tion of the WD method to larger and spatially more extended 
surveys (RAVE, Gaia). 

5. Age-kinematics characterization 

Two different approaches are adopted in this section. First, the 
whole sample is divided by age range into statistically signif- 
icant subsamples to trace the kinematic structures in the U—V 
plane through age. Secondly, only those stars with the most 
precise ages -relative errors of less than 30%- are used to 
study structures and periodicities in U-V-age space and other 
connections between kinematic structures and the evolutionary 
state of their members in the context of the branches proposed 
in Sect. 14.21 Note that in both analyses the sample of stars used 
only includes the FGK and OBA dwarfs samples, since individ- 
ual ages are not available for the other subsamples (see Table 

5.1. Age dependence in the U-V plane 

Figure Q~3] shows the distribution in the U—V plane and the dis- 
tribution for the V rot component obtained by WD for subsam- 
ples with different age ranges. The age bins have been cho- 

17 In Sect. [2] we mention that the Ho- equivalent width could be used 
as information on age for the M dwarfs. We have checked, however, 
that the number of stars per sample when active and non-active stars 
are analysed separately is too low to derive significant results from 
the distributions: too smoothed distributions are obtained by WD with 
J pl ateau- This study has been postponed until the use of recent new data, 
i.e. lBochanski et alj d2007» . 
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Fig. 13. Left column: Density field in the U-V plane obtained 
by WD with J - A for subsamples of different ages (in Gyr): 0- 
0.1 (1792 stars), 0.1-0.5 (1501), 0.5-2.0 (3368), 2.0-4.0 (3917), 
4.0-8.0 (2561) and >8.0 (2053). Right column: Density field 
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V^kms" 1 ) oge (Gyr) 

Fig. 14. Left: Density field in the V rot -age plane of the stars with e age < 30% (7016 stars) obtained by WD. Right: Integrated 
density of stars in each branch worked out using the WD method as a function of age for stars with e age < 30% (223, 1 195, 1078 
and 806 stars in the Hercules, Hyades-Pleiades, Coma Berenices and Sirius branches respectively). 



sen in order to emphasize the changes in these distributions. 
Important observations derived from Fig. [13] are: 

- For the youngest stars, only the Pleiades and Coma 
Berenices kinematic groups appear and the branches are 
not identified in the U—V plane. For ages between 0. 1 and 
0.5 Gyr, the Hyades and Sirius structures have begun to ap- 
pear. The Hyades structure is the most prominent in all the 
subsamples with age > 0.5 Gyr. 

- Certain structures appear and disappear (or at least their 
contrasts vary considerably in relation to other structures) 
when subsamples with different age ranges are consid- 
ered. See, for instance, the Pleiades structure at (U, V) ~ 
(-12, -22)kms~ 1 . From top to bottom, we see how this 
structure is clearly evident for very young stars but is less 
significant in the 0.5-2.0 Gyr range, only to reemerge at 4-8 
Gy0 

- The structure of branches is well traced up to the oldest 
samples. The minimum age of the stars in each branch is 
discussed in Sect. 15.21 

- The separation between branches does not depend on age. 
The value of 15kms _1 is found to be appropriate except 
for the youngest stars (<100Myr) where the Pleiades and 
Coma Berenices moving groups are closer in V rot . 

- The drop of density at the left of the Hyades-Pleiades 
branch considered in Sect. I4.2l is noticed in all age ranges. 
In addition, we observe that the relative density between 
this branch and the Hercules branch clearly decreases with 
age. 



18 lAsiain et alj dl999bl) . working only with the sample of B and A 
main sequence stars, found a significant variation of the Pleiades mov- 
ing group with age. However, as ages in their sample are < 0.5 Gyr, 
their interpretations cannot be applied here. 



5.2. Structures and periodicities in the V rot -age plane 

Figure [14] (left) shows the distribution in the V ro ,-age plane 
for the stars with well-defined ages (e age < 30%, 7016 stars). 
The age distribution of each branch (Fig. [14] right) is com- 
puted from the whole distribution, selecting the region of the 
branches as detailed in Sect. l4.2l (centres at -55, -25, -10 and 
5 kms _I for the Hercules, Hyades-Pleiades, Coma Berenices 
and Sirius branches respectively with a width of +4kms~ 1 ). 
Important observations derived from this distributions are: 

- An extended age distribution is confirmed for all four 
branches. 

- While the Hyades-Pleiades and Coma Berenices branches 
have an important fraction of very young stars, Sirius has 
its first main peak for slightly older stars at ~ 400 Myr. 
Consequently, the velocity distribution becomes structured 
by the tree mai n branches for stars of ~ 400 Myr. 



Dehnenl (120001) assigns an age of > 8 Gyr to the Hercules 



group and iBobylev & Bajkoval (12007) claim that it does 
not show any branch-like shape unless very old (> 8 Gyr) 
subsamples are studied. Here the kinematic region of the 
Hercules branch is populated with stars of ages of more 
than 1 Gyr and its extended branch-like shape is seen in all 
subsamples with ages > 2 Gyr (Fig.[T3l> where it reaches to 
a quite extended longitude. 

A clumpy distribution is observed inside the Hyades- 
Pleiades branch, with a periodicity in age of about 500-600 
Myr, maintained at least in the 0-2.5 Gyr range. However, 
this periodicity is marginally significant as the abso- 
lute error in age of some stars exceeds the period ob- 
tained. We point out that this periodicity should be con- 



19 Notice that the periodicity is statistically significant in the sense 
of the discussion in Sect[3](-/ = J plateau)- 
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tr asted in the future with the star formation rate obtained 
in Hernandez et al. (2000), where an oscillatory component 



of period ~ 500 Myr is found. For the other branches only 
an outline of the shape of the whole age distribution is ob- 
served. 



6. An attempt to evaluate [Fe/H] dependence 

The study of the [Fe/H] metallicity distribution of the stars 
along the branches and its relation to kinematics and age is re- 
stricted to the FGK dwarfs, which is the only sample for which 
metallicity data is available (see TableQ]). In Table|4]we present 
the mean metallicity and dispersion values for each branch with 
stars selected as di scussed in Sect. 14.21 The values obtained by 
Helmi et alj (|2006), using the same sample but a different def- 
inition of the position of the kinematic structures, are included 
in the table for comparison (Coma Berenices was not included 
in Helmi's analysis). These authors found a common metallic- 
ity dispersion of 0.2 dex for all three superclusters. However, 
we obtain slightly lower mean metallicities for all the branches 
and a higher m etallicity dispersion for the Hercules branch. 
Hayw ood (2006) point ed out that systematic biases in the data 
from INor dstrom et al. (2004) provoke an artificial increase in 
the [Fe/H] dispersion at least for ages < 3Gyr. As the four 
branches present a wide range of ages (Sect. [5]), it will be as- 
sumed that this specific bias influences each of the branches 
similarly. 

More importantly, in Fig. [15] we present the distribution in 
the [Fe/H]-V rot plane. Both the WD method and the rotation 
of the velocity components by 16° allow us to identify the four 
conspicuous branches and their rather wide range of metallic- 
ity. A general correlation can be seen: a more negative V rot im- 
plies a higher mean metallicit\F'l. Notice that the correlation is 
complex due to either the bimodality of the Hyades-Pleiades 
branch or to the fact that the Hercules branch does not follow 
the overall pattern of the three main branches. 

The bimodality for the Hyades-Pleiades branch shows the 
peaks at [Fe/H] ~ -0.05 and [Fe/H] ~ -0.15. The metallic- 
ity distribution for the stars along the position that this branch 
traces in the U—V plane is studied in detail in Fig. [16] We ob- 
serve that the [Fe/H] distribution evolves from the Hyades 

{Vrot SOkms 1 ) with the peak at [Fe/H] ~ -0.05 - 

in agreement with the results of iFamaey et al.l (120071) - to the 
Pleiades, which placed at U rot ~ -5kms _1 shows a peak at 
[Fe/H] ~ -0.15. 



Figure [17] (top) shows the distribution of the whole sample 
in the age-[Fe/H] plane. Due to the complexity in the age- 
met allicity relation (see the exhaustive discussion on the biases 



in 



Havwoodl 2006) any new conclusion is beyond the scope 



of this paper. Our analysis is restricted to a comparison be- 
tween branches. Figure [TT] shows the distributions in the age- 
[Fe/H] plane for the stars along each of the four branches. 
To a first approximation these distributions exhibit the same 



20 Notice that this c orrelation was not as clearly observed in 
INordstrom et alj J2004) (see their Fig. 32), despite using the same 
sample, which demonstrates the capabilities of the WD method and 
the rotation. 
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Fig. 15. Density field in the V mt -[Fe/H] plane for the sample 
with available [Fe/H] (13109 stars) obtained by WD. 

0.000 0.070 0.141 0.211 0.281 0.352 0.422 0.493 
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Fig. 16. Density field in the U rot -[Fe/H] plane obtained by 
WD for the Hyades-Pleiades (2271 stars). 

Table 4. Mean metallicities and dispersions for the branches 
considered in the present work and for the superclusters ac- 
cording to (l) lHelmi et~ai1(l2006l) . 





Branches (this work) 


Superclusters (1) 




[Fe/H] <r lFe / H ] 


[Fe/H] o- [Fe/ n] 




( dex) ( dex) 


( dex) ( dex) 


Hercules 


-0.15 0.27 


-0.13 0.2 


Hyades-Pleiades 


-0.11 0.20 


-0.08 0.2 


Coma Berenices 


-0.16 0.22 




Sirius 


-0.21 0.21 


-0.18 0.2 



general tendency as the whole sample. The WD treatment al- 
lows us to detect other significant features such as the clumps 
in age for the Hyades-Pleiades branch detected in Sect. [5] which 
show a decrease in [Fe/H] with increasing age. For the range 
1-3 Gyr, where more data is available, although it is prelimi- 
nary, the slope of the Sirius distribution seems to be slightly 
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0.0043 0.0050 
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age (Gyr) 

0.027 0.036 0.045 0.054 



0.022 0.029 0.036 0.043 0.050 



Fig. 17. Density field in the age-[Fe/H] plane for all stars with 
available metallicities and ages (11215 stars) and for the differ- 
ent branches of Hercules, Hyades-Pleiades, Coma Berenices 
and Sirius (436, 1973, 1559, 1229 stars respectively) obtained 
by WD. 



more negative than for the other two main branches (from 
~ -0.10 ± OMdex/ Gyr for Sirius to ~ -0.07 + 0.01<iex/ Gyr) 
for Hyades-Pleiades. 



7. Summary and discussion 

Moving groups are becoming a powerful tool for studying the 
large-scale structure and dynamics of the Milky Way. We ap- 
ply multi-scale techniques -wavelet denoising- to an extensive 
compendium of more than 24000 stars in the solar neighbour- 
hood to characterize the observed kinematic structures in U—V— 
age-[Fe/H] 4-dimensional space. In this paper we focus our 
analysis on establishing the observational constraints that will 
allow us to study the origin and evolution of these structures. 

The advent of Hipparcos astrometric data led to the defini- 
tive recognition of a non-smooth distribution function of the 
Galactic disc. Our results corroborate this and go one step fur- 
ther towards characterizing the velocity distribution function. 
Branches connecting the classic moving groups in the U—V 
plane are definitely the dominant structures. We confirm the 
existence of the Sirius, Coma Berenices, Hyades-Pleiades and 
Hercules branches. The first three branches are spaced at in- 
tervals of approximately 15 kms -1 with no significant varia- 
tions with age or spectral type. The Hercules branch is located 
about 30 kms -1 from the Hyades-Pleiades branch. The four 
branches present a negat ive slope of /3 ~ 16° in the U—V plane, 
lower than that found by Skuli an et al. (1999). We have studied 
the density drops in the U—V plane and the existence of abrupt 
edge lines is corroborated, especially at low U and high V and 
between the Hyades-Pleiades and the Hercules branches. Each 
branch presents a different density distribution in its extremes 
near these edge lines. These features definitely rule out the clas- 
sic idea of a smooth velocity field distribution. The use of the 
same statistical method for the samples with different spec- 
tral types has made possible an exhaustive comparison between 
their kinematic planes. The branches are present in all the dis- 
tributions. The study of variations of the kinematic structures 
with Galactic position has been presented for the KM giant 
sample. A significant change of contrast among substructures 
inside the branches is confirmed and the shape of the Hercules 
branch changes among regions, being more conspicuous in the 
region towards the anti-rotation direction with respect to the 
centre of this sample. 

We have analysed, for the first time, the age and metal- 
licity distributions of the branches. We confirm an extended 
age distribution for all the branches but a different mini- 
mum age of their stars. The set of the three main branches 
of Hyades-Pleiades, Coma Berenices and Sirius is well es- 
tablished for stars > 400 Myr and the extended branch-like 
shape of Hercules is detected in all subsamples with ages 
> 2 Gyr. The relative density between the Hyades-Pleiades and 
the Hercules branches clearly decreases with age. We find a pe- 
riodicity in age of about 500-600 Myr in the Hyades-Pleiades 
branch. For the other branches only an outline of the shape of 
the whole age distribution is observed. A wide range of metal- 
licity is found for each branch, especially for Hercules with 
a higher metallicity dispersion. A complex relation between 
kinematics and metallicity has been established. Concerning 
the age-metallicity relation, differences are observed between 
branches. The periodic bumps in the age distribution for the 
Hyades-Pleiades branch show a decrease in metallicity with in- 
creasing age. 
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All the above observational results and, moreover, the 
whole set of distributions in the U-V-age-[F e/H] space pre- 
sented here, are the fundamental elements necessary to check 
the present and future dynamic models proposed for the for- 
mation of kinematic structures. While some of these cur- 
rent models already explain some of the observational results, 
other important observational f eatures remain unexpla ined. 
For instance, the simulations of iDe Simone et aD d2004l) pro- 
duce branches with a slop e which fits that measure d here. 
Furthermore, the model in Oui llen & Minchev (2005) agrees 
with the fact that the Hyades-Pleiades branch is more metal- 
lic than the Coma Berenices branch. However, features such 
as the mix of stars in the Hyades and Pleiades moving groups 
in periodic clumps in age in the same branch -perhaps reflect- 
ing a common mechanism for the formation of these two kine- 
matic structures- or the peculiar metallicity distribution inside 
branches demand an in-depth explanation. 

To make further progress in this field, simulations obtained 
from orbit integration under a model for the Galactic potential 
are being undertaken (Antoja et al. 2008, in preparation). After 
checking that the simulated kinematic structures critically de- 
pend on the gradients of the spiral arm potential, two differ- 
ent models for the Galactic spiral arms are being explored: the 
classic cosine perturbation ( Linl 1 197 II) and a 3D mass distri- 
buti on (superposition o f the inhomogeneous oblate spheroids 



Pichardo et al. 2003). We are also incorporating stellar ages 



and metallicities in the simulations to trace the origin and evo- 
lutionary state of the test particles, which is another key ele- 
ment. Of course, other ingredients such as the Galactic bar will 
also be included to analyse the kinematic structures in an over- 
all dynamic context. 

Thus, the present detection and characterization of the 
branches in U-V-age-[Fe / H] 4-dimensional space provide us 
with the main properties summarized above. It is the whole set 
of distributions in this space that contains all the information. 
The application of the same statistical technique to both the 
observed and simulated data will allow a direct comparison be- 
tween them, and thus be a powerful test of the models for the 
formation of the structures. This comparison will eventually of- 
fer a physical interpretation of the formation of the kinematic 
structures in terms of local or large-scale dynamics. 
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